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Abstract. For many complex systems the interaction of different scales is among 
the most interesting and challenging features. It seems not very successful to extract 
the physical properties in different scale regimes by the existing approaches, such as 
structure-function and Fourier spectrum method. Fundamentally these methods have 
their respective limitations, for instance scale mixing, i.e. the so-called infrared and 
ultraviolet effects. To make improvement in this regard, a new method, multi-level 
segment analysis (MSA) based on the local extrema statistics, has been developed. 
Benchmark (fractional Brownian motion) verifications and the important case tests 
(Lagrangian and two-dimensional turbulence) show that MSA can successfully reveal 
different scaling regimes, which has been remaining quite controversial in turbulence 
research. In general the MSA method proposed here can be applied to different 
dynamic systems in which the concepts of multiscaling and multifractal are relevant. 

Keywords'. intermittency, multifractal, two-dimensional turbulence, Lagrangian 
turbulence 


Multi-level segment analysis: definition and application in turbulent systems 

1 . Introduction 


2 


Multiscale is one of the most important features commonly existing in complex systems, 
where a large range of spatial/temporal scales coexist and interact with each other. 
Typically such interaction generates scaling relations in the respective scale ranges. 
To understand the multiscale statistical behav ior has been remaining as the research 
focus in various areas, such as fluid turbulence flFrischl. Il995l). hnancial market analysis 


flMantegna and Stanley 


science fISchmitt et ah 


19961: ISchmitt et al.l. Il999l: iLi and Hnand. 12014 ) , envi ronmental 


19981) and population dynamics fISeuront et ah 


19991) . to list a 


few. Among a number of existing analysis approaches for such problems, the standard 
and most cited one is structure-function (SF), which is hrst introduced b y Kolmogorov 
in hi s famou s homogeneous and isotropic turbulence theory in 1941 (jKolmogorov . 


194ll: iFrischl . 119951) . However, the average opera tion in SF mixes regions with 
different correlations flWang and Petersl . l2006l. 120081) . Mathematically, SF acts as 
a hlter with a weight function of W(k, £) = 1 — cos(2nki), in which k is the 


wave number and £ is the separation scale flDavidson and Pearsonl . l2005l: iHuang et ah 


20101 ). It thus makes the statistics at different scale £ strongly mixed, resulting in 
the so-called infrare d and ultraviolet e ffects, respectively for large-scale and small- 
scale contamination flHuang et al.l, 120131) . The situation will be more ser ious when an 


energ etic structure presents, e.g., annual cycle in collected geoscience data flHuang et ah 


20091 ). large-scale circul ations in Ray l eigh-Benard convec tion, vortex trapping events in 


Lagrangian turbulence flHuang et al.l . 120131: IWangl . 120141) . 

Fourier analysis in the frequency domain has the similar dehciency as SF, i.e. any 
local event will propagate the influence over the entire a nalyzed doma i n, especiall y 
for th e nonl inear and nonstationary turbulent structures (IHuang et al.l Il998l . 120111 ). 


Fargd fjl992l) claimed that a localized expansion should be preferred over unbounded 


trigonometric functions used in Fourier analysis, because it is believed that trigonometric 
functions are at risk of misinterpreting the characters of held phenomena. An 
alternative approach, namely wavelet transform is then proposed to overcome the 


possible shortcorn ing of the Fourier transform with local capability (jPaubechiesl. Il992 


Farge et al.l. 119961) . However, the same problem as Fourier analysis still exists, if the 


hxed mother wavelet has a shape different from the analyzed data structure. We also 
note th at the classical structure-fun ction analysis is referred to as ‘the poor man’s 
wavelet’ (jLovejoy and Schertzerl . 120121 ) . 

To overcome the potential weaknesses of SF or Fourier analysis, several 
methodologies have been proposed in recent years to emphasize t he local geometrical 


features, such a s wave let-based methodologies (wavelet leader fl.Taffard et al. 

wavelet transform modulus maxima f Muzv et al. 


Lashermes et al. 


2008’ 


Oswiecimka et ah. 

2001 

3), 


detrended structure-func tion (iHuand. 120141). scaling of maxim um probability density 


analysis flHuang et al. 


20081 


Iboiil), 


huct nation analysis (jPeng et ahl . Il9941 ). 


2005 


1993 


function of increment s JHuanget ah . 2011 : Qiu et ah . 2014 ). and Hilbert spectral 


to name a few. Note that different approaches 
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may have different performances, and their own advantages and disadvantages. For 
example, the detrended structure-function can constrain the influence of the large- 
scale structure, using the detrending procedure to remove the scales larger than the 
separation distance i. I n practi c e, the famous 4/5-law can then be more clearly retrieved 
than the classical SF flHnangi . I2ni4l ). However, this method is still biased with the 
vortex trapping event in Lagrangian tur bulence, which typi cally possesses a time scale 
around 3 ~ 5r^ in the dissipative range flToschi et al.l. 120051) . The scaling of maximum 
probability density function of increments helps to quantify the background fluctuation 
of turbul ent fields. Compared with SF , it ca n efficiently extract the first-order scaling 
relations flHnang et al.l. 1201 luQin et al.l. 120141 ): however, it is difficult to extend to higher 
gth-order cases. 

A new view on the held structure is based on the topological features of the 
extremal points. In principle, physical systems may assume different complex i ty ari d 
interpretability in different spaces, such as physical or Fourier (jWang and Petersl . 120131) . 
The extremal point structure in physical space has the straightforward advantage 
in dehning characteristic parameters. Considering a huctuating quantity, turbulence 
disturbs the how held to generate the local extrema, while viscous dihusion will 
smooth the held to annihilate the extremal points. By nature the statistics of local 
extremal points inherits the process physics. Based on this idea, Wang, Peters and 
other co llaborators have studied pas si ve sca l ar tu rbulence via dissipation element 


analysis flWang and Petersl . 


12006 


20081) . IWangl (120141) analyzed the Lagrangian velocity 


by dehning the trajectory segment structure from the extrema of particles’ local 
acceleration. Such diagnosis verihes successfully the Kolmogorov scali ng relation, which 
has been argued controversially for a long time (IFalkovich et al.l . 120121 ). However, under 
some circumstances the extremal points may largely be contaminated by noise, thus 
partly be spurious. In other words extremal points are sensitive to noise perturbation. 
Although data smoothing can relieve this problem, some artihcial arbitrariness will 
inevitably be introduced; moreover it may not be easy to design reliable smoothing 
algorithms from case to case. In this regard the extrema-based analyses are not generally 
applicable, e.g. with noises from measurement inaccuracy, interpolation error or external 
perturbations. 

In this paper a new method, multi-level segment analysis (MSA), has been 
developed. The key idea hereof is based on the observation that local extrema are 
conditionally valid, indicating a kind of multi-level structure. Compared with the 
aforementioned extrema-based analyses, this new method is a reasonable extension with 
more applicability. Details in algorithm dehnition, verihcation and applications will be 
introduced in the following. 


2. Multi-level segment analysis: method definition 

Considering any function f{x) in some physical process, where x is the independent 
variable, e.g., the spatial or temporal coordinates, xq is a local extremal point with 
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s=100 



Figure 1. (Color online) An illustration of the extracted extreme points at different 
s levels: a) s = 100, and b) s = 200. The local maxima and minima are demonstrated 
respectively by O and □. The vertical line indicates the window size s. To ensure 
spacial homogeneity, a sliding window with level s scans over the whole data set, see 
more detail in the text. 


respect to scale s is defined as /(xq) < f{x),\/x G (xq — s,a:o + s) (minimum), or 
f{xo) > f{x),\/x G (xo — s,xo + s) (maximum). Extrema are conditionally valid. For 
instance, if Xq is extremal at scale s, it may not be extremal at a larger scale Si > s. 
Figure [U illustrates that for an artificially generated signal, at different s levels both 
the number density and the fluctuation amplitude of the extremal points will change 
accordingly. Under some special conditions extremal points have simple but interesting 
properties. For instance, for a monotonous function, there is no extrema for any s; for a 
single harmonic wave function with a period of T, the number of extremal point remains 
constant. For real complex multiscale systems as turbulence, variation of extremal point 
is continuously dependant on s. 

At a specihc s level, denote the corresponding extremal point set as i = 1, 2,... 
(along the coordinate increasing direction). Numerical tests show that typically these 
points are max — min alternated for small s, while when s increases max — max or 
min —min events may also appear but with low probability (e.g., 5%), depending on 
the /(x) structure and the process physics. The segment is dehned as the part of /(x) 
between two adjacent extremal points. The characteristic parameters to describe the 
structure skeleton are the function difference, i.e., /(xs,i) — f{xs^i-i) and the length 
scale, i.e., E = Xs^i — Xs,i-i. By varying the s value, different extremal point sets, and 
thus different segment sets, can be obtained. In this sense this procedure is named 
as multi-level se gment analysis (MSA); while the existing approaches based on local 


extremal points flWang and Petersl . 120061120081: IWangi. 120141 ) calculate extrema from the 
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DNS data; thus can be understood as single-level. Compared with the conventional SF, 
in MSA the segment length scale is not an independent input, but determined by f{x). 

For a specified s, scan over the data domain to have the corresponding segment 
characteristics, i.e. f{xs,i) — /(xs^j_i) and i = Xs,i — Xg^i-i- Collect the results for 
different s, the functional statistics can be defined. Numerically the same segments 
may be counted repeatedly for different s, which then need to be excluded. In terms of 
structure function, the mathematical expression is (for the gth order case) 


- f{Xs,i-l)]'^\Xs,i - Xs,i-1 = i)s, ( 1 ) 

where {■)s denotes sampling over different s. If any scaling relation exists, one may 
expect a power-law behavior as 

Dg(£)(2) 

in which ^{q) is the MSA scaling exponent. 

For comparison we also include here briefly the classical SF and wavelet-leader 
definitions. The gth order SF is written as 

Sfii) = {\Axe{m^f^^\ (3) 

where Axfit) = x{t + i) — x{t) and £ is the time separation scale. The scaling exponent 
Q{q) characterizes the fluctuation statistics, (^{q) is linear for monofractal processes such 
as fractional Brownia n motion, and nonlinear and concave for multifractal processes 
(ISchertzer et al.l . 119971 ). As mentioned above, SF mixes information from different scales. 
It is al s o lirn i ted by the slope of the Fourier spectrum E{f) ~ f~^, e.g., 1 < /3 < 3 
flFrischl . Il995l: iHnang et al.l . 120101 ). We denote this as /5-limitation. 

for example, w a velet- 


There are several different wavelet-based methoc 


transform-inodulu s-maxima (WTMM) llMnzv et ah 


Lashermes et ah 


Muzv et ah . 1993) and wavelet-leader (WL) f Jaffard et ah . 20051: 


1991 


s, 


Mallat and Hwang 


Wendt et ah 


1992 


20071: 


200811 . We consider here only WL. More detailed discu s sions of those 


meth ods can refer to Ref. fjjaffard et al.l. l2005l: lOswiecimka et al.l. 120061: iHnang et al. 
201111 and references therein. 

The discrete wavelet transform is defined as 


'fi{k,j)= / x{t)(p{2 H-k) df. 


( 4 ) 


where ip is the chosen wavelet, 'fi{k,j) is the wavelet coefficie nt, k is the position index 


?' is t he scale index, and E = 2^ is the corresponding scale (IDaubechiesl . Il992l: iMallat 
19991 1. Every discrete wavelet coefficient 'il!{k,j) can be associated with the dyadic 
interval g{k,j) 

g{kj) = [2^k,2^{k + l)). (5) 

Thus the wavelet coefficients can be represented as = 'ip{k,j). The wavelet-leader 
is defined as 


l{kj)= sup li^ig')], 

g'c3g(k,j),j'<j 


( 6 ) 
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where 3p(k, j) = pjk — 1, j) U Q{k, j) U g{k + 1, j) (jjaffard et al.l . l2005l: iLasliermes et all 


2008 


Wendt et all 120071) . The expected scaling behavior can be expressed as 


(7) 


in which x(g) is the correspon ding scaling exponent. The calcnlation efficiency has been 


discu ssed for various datasets (IJaffard et ah 

2oS). 


20051: 1 Wendt et ahl. 120071: iLashermes et ah 


As shown in the rest of this paper, for simple cases, MSA and the classic methods 
show pretty identical results; while for complex analyses as turbulence, because of the 
algorithmic principle to depict the function structure, the former one is more effective 
and efficient. 

Some additional comments are stated as follows. First, MSA can be considered as 
a dynamic-based approach without any basis assumption a priori. Since local extrema 
imply the change of the sign of f\x) across Xs^i- Assuming that f{x) represents a velocity 
signal, f'{x) is then the acceleration, i.e. a dynamical varia ble. Second, MSA shares the 
spirit of the wavelet-transform-modulus-maxima (WTMM) flMnzy et al.l . ll993h . in which 
only the maxim um modulus of the wavele t coefficient is considered. However, as argued 
in several Refs. fiHuang et ahl. Il998l. l201ll) . if the shape of the chosen mother wavelet is 
different with the specihc turbulent structure, e.g., ramp-cliff in the passive scalar held, 
additional high-order harmonic components are then mixed to ht the difference between 
the physical s t ructu re and the mother wavelets, which then biases the extracted scaling 


flHuang et al.l . 1201 ll) . In this aspect MSA can be considered as a data-driven type of 


WTM M without any transform. Similar with the wavelet leader, IWelter and Esquel 
(120131 ) proposed a multifractal analysis based on the amplitude extrema of intrinsic 
mode functions, whi ch can be retrieved by the empirical mode decomposition algorithm 


flHnang et al.l . Il998l ). In some synthesized data tests, a proposed parameter m = 2 is 


used to determine the search domain for the local amplitude maxima. 


3. Case verification and applications 


3.1. Fractional Brownian motion 


Fractional Brownian i notion ffBm) is a ge neralization of the classical Brownian motion. 
It was introduced by iKolmogorovI (119401) and extensively studied by Mandelbrot and 
co-workers in the 1960s (IMandelbrot and Van Nessl . |l968|) .Since then, i t is con si dered 
as a classical scaling stochastic process in many fields (IBeranl . 1 1994 iRogersl . Il997l: 


Donkhan et ahl. 120031 ). In the multifractal context, fBm is a simple self-similar process. 


More precisely, the measured SF scaling exponent Ch{q) is linear with the moment 
order q, i.e., Ch{<i) = qH, in which H G (0,1) is the Hurst number. The above 
linear scaling relation has been verihed by various methodologies, such as, classical SF, 
wavelet-based method, detrended fluctuation analysis and detrended structure-function, 
to list a few. In the present work, a fast Fourier transform based Wood-Chan algorithm 
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Figure 2. (Color online) a) Calculated second-order I? 2 (^) by MSA for fractional 
Brownian motion with Hurst number H = 0.2 (o), H = 1/3 (□), H = 0.5 (v) 
and H = 0.75 (A), respectively. Power-law behavior can be observed for all H. b) 
Compensated curves 772(£) x . For display clarity, these curves have been vertical 
shifted, c) Experimental scaling exponents (?) in tbe range 0 < g < 4. The errorbar 
implies the standard deviation from 100 realizations, d) Measured singularity spectrum 

/(a)- 


(IWood and Chad . Il994l) is used to synthesize the fBm data with 100 realizations, and 
10® points for each H. 

Figure |2]a) shows respectively the calculated second-order 'D 2 {t) by MSA for Hurst 
number H = 0.2 (o), H = 1/3 (□), H = 0.5 (v), and H = 0.75 (A). The power-law 
behavior is observed for all H considered here. To emphasize the agreement between the 
measured and the theoretical value, compensated curves of are shown 

in hgure [2]b). For display clarity, these curves have been vertical shifted. Visually, when 
H < 1/3, the measured ^{q) deviates from the theoretical prediction. Such outcome may 
be related to MSA itself or the fBm date generation algorithm. Further investigation 
will be conducted hereof. Figure |2]c) shows the measured ^{q), which are estimated in 
the range 10 < £ < 10, 000 by least-square-htting. The errorbars indicate the standard 
deviation from 100 realizations (same in the following). These curves demonstrate for 
all the cases the linear dependence of the measured scaling exponent ^{q) with q, whose 
slope is H. In other words, multifractality can successfully be detected by MSA for all 
H values (including H < 1/3). Consider the so-called singularity spectrum, which is 
dehned through a Legendre transform as, 

d^(?) 


a = —/(«) = min - ^(g) + 1} • 

dq q 


( 8 ) 


For a monofractal process, a is independent with q, e.g., a = H, and /(a) = 1 (IFrischl . 
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Figure 3. (Color online) Experimental results for the SF analysis: a) the measured 
second-order SFs <S'2(t); b) compensated curves 5*2 (r) x c) measured scaling 

exponent C(9); and d) the measured singularity spectrum /(a) versus a. Symbols are 
the same as in figure [2j 



Figure 4. (Color online) Experimental results for the wavelet leader: a) the measured 
second-order Z 2 {£)', b) compensated curves Z 2 {£) x c) measured scaling 

exponent xia)i and d) the measured singularity spectrum /(a) versus a. The symbols 
are the same as in figure [5] 


1995h . In practice, for a prescribed q, a wider range of a has, a more intermittent the 
process is. Figure |2]d) shows the measured /(a) in the range 0 < g < 4. It conhrms the 
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Figure 5. (Color online) Comparison of the measured Hurst number H versus the 
given H. The theoretical value is indicated by the the solid line. 


monofractal property of the fBm process. 

Figure [3] lists the results from SF: a) the measured second-order SFs, 5*2 (r), b) the 
compensated curves 5'2(r) x , c) the measured scaling exponents ({q) and d) the 
corresponding singularity spectrum /(a) versus a, respectively. Visually, a plateau is 
observed for all r, showing a perfect agreement between the detected scaling and theory. 
Here the scaling exponents ({q) are also estimated in the same the range 10 < r < 10, 000 
by least-square-fitting. It shows that C{q) is linear with q, and the measured singularity 
spectrum f{a) detects correctly the monofractal property of the fBm process. 

Figure H] shows the results from the second-order WLs. Power-law behavior is 
observed for all H. Compensated curves Z 2 {i) x £-(2^^+2) gPow a clear plateau when 
^ > 10. Similar as MSA, misalignment is observed when ^ < 10, which may be due 
to the fBm generator used in this study. Scaling exponents are then estimated in 
10 < r < 10,000. The measured x(g) and the corresponding singularity spectrum /(a) 
capture the monof ractal proper t y of the fBm. Note that the singularity spectrum of 
WLs is defined as flHuang et al.l . l201ll ) , 


a = 


dg 


- 1. /(o) = mill {o? - x(«) + 1 + 9} 


( 9 ) 


Figure [5] shows the comparison of the different estimated Hurst numbers H, which 
are calculated by linear htting of the measured scaling exponent i{q)^ C(?) s-iid x(g), 
respectively. Visually, SF and WL provide almost the same performance, while MSA 
slightly overestimates H when H > 1/3. All methods considered here confirm the 
monofractal property of the fBm process. Considering th at the turbule n t data are much 
different from the simple fBm, as already shown in Ref. iHuang et al.l (120111) . both SF 
and WL are strongly influenced by the real turbulent structures. 

We provide a comment here on the deviation of the measured ^(g) from the 
theoretical prediction when H < 1/3. In the context of extreme point based MSA, the 
intrinsic structure of fBm is presented by these extrema. Therefore how ^{q) will behave 
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is determined by the relation between the Hurst number H and the distribution of the 
extreme points, as well as the fBm data generation. Or in other words, the retrieved 
^(g) relies on the dynamic behavior of the process itself, which is deeply related with 
the distribution of the extrema point. This is beyond the topic of this paper. We will 
present more details on this topic in the follow-up studies. 

3.2. Lagrangian velocity 



10^ 


10 “ 


10 “ 


I I I mill—I I I mill—rm 


it™, 


,V VVVVv, 


b) 


- cfjnniiiiiujjV 

—- V -^ 

V □ A o 

cv ° o 


L A 


.I . .I _l_LU 


0.1 


10 100 1000 


Figure 6. (Color online) a) Results from various methods: classical SF (v) and 
the Hilbert-based method (o) for the second-order statistics, the trajectory-segment 
method (A) and MSA II’i(t) (□) for the first-order statistics. For the Hilbert-based 
statistics, frequency has been converted to time by t = l/w. b) The corresponding 
curves compensated by the dimensional scalings, i.e. S 2 {t) ~ £t and 5 'i(t) ~ 

There is no plateau from SF, in consistency with other reports in the literature. The 
following convincing scaling ranges can be observed: about 2 < r/r^ < 60 from MSA, 
2 < r/r,, < 50 from the trajectory segment analysis and 10 < r/r,, < 100 from 
the Hilbert-based method, which verify the prediction of the Kolmogorov-Landau’s 
phenomenological theory. For display clarity, these curves have been vertical shifted. 


The Lagrangian velocity SF has been extensively studied. Because the time scale 
separation in the Lagrangian frame is more Reynolds number dependent than the length 
scale case in the Eulerian frame, the hnite Reynolds number influence b ecomes stronger. 


maki ng the Lagrangian velocity scaling relation quite controversial (iFalkovich et ah 


20121 ). More specihcally, this is recognized as a consequence of mixing between 


l arge-scale structures and energet i c sma ll-scale structures, e.g., vortex trapping events 
(iToschi et al.l . 120051: iHnang et al.l . l2013l ). The gth-order SF of the velocity component 
Mj (i = 1, 2, or 3) is dehned as: 

Sfir) = {[ufit + r) -Ui{t)fi), (10) 

where r is an arbitrary tir ne separation scale. Fr om dimensional analysis, the 2nd-order 
SF is supposed to satisfy (IFalkovich et al.l. 120121) : 

S2{t) = Coer, ( 11 ) 


where e is the rate of energy dissipation per unit mass and Cq is assumed as a universal 
constant at high Reynolds numbers. 
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To analyze this problem, we adopted the data from direct numerical simulation 
(DNS), implemented for the isotropic turbulence in a 2048^ cubic domain. The boundary 
conditions are periodic along each spatial direction and kinetic energy is continuously 
provided at few lowest wave number components. A hne resolution oi dx ^ rj (the 
Kolmogorov scale) ensures resolving the detailed small-scale velocity dynamics. The 
Taylor scale A based Reynolds number Re\ is about 400. Totally 0.2 million Lagrangian 
particle samples are collected, each of which has about one integral time life span. 
During the evolution process, the velocity and velocity derivatives are recorded at 
each r „/20, which 7 y is t he Kolmogorov time. More numerical details can be found 


in Ref. iRenzi et a 1.1 (120091 ) and references therein 


Re cently , this database has been analyzed respectively bv lHuang et al.l (120131) . and 


Wangl (I 2 OI 4 J ) to identify the inertial scaling behavior. The former study employed the 
Hilbert-based approach, in which different scale events are separated by the empirical 
mode decomposition without any a priori basis assumption and the corresponding 
frequency u is extracted by the Hilbert spectral analysis. They observed clearly an 
inertial range of 0.01 < < 0.1, i.e. 10 < r/r,, < 100. The scaling e x ponen ts C(?) 

agree well with the multifractal model (see details in Ref. iHuang et al.l (120131 )). The 
latter one studies the extrema of the fluid particle acceleration, which physically can be 
considered as the boundary markers between different flow regions. With the help of 
the so-called trajectory segment structure, the clear scaling range does appear. Because 
of int erpolation inaccuracy (noise), DNS data need to be particularly smoothed (IWand. 
2 OI 4 J ). which may lead to some artihcial input. 

Figure [6]a) shows the numerical results from various methods: classical SF (y) and 
the Hilbert-based method (O) for the second-order statistics, the trajectory-segment 
method (A) and MSA (□) for the hrst-order statistics. Except for SF, clear power-law 
behaviors can be observed for others. To emphasize this, curves compensated by the 
dimensional scalings, i.e. S 2 {t) ~ er and 5 'i(t) are plotted in Fig.|6]b). The 

SF curve does not show any plateau, which is co nsistent with reports in the literature 
( Falkovich et ah . 2012 : ISawford and Yenii^ . 2011 ). even for high Re cases as Re\ ~ 815 
experimentally or Rex — 1000 numerically. The absence of the clear inertial range makes 
the Kolmogorov-Landau’s phenomenological theory quite controversial. In comparison, 
the MSA curve shows a convincing plateau in the range 2 < r/y < 60. 

It need to mention that for this DNS database the inertial range has been recognized 
as 10 < t/t„ < 100 for both t he single particle statistics using the Hilbert-based 


methodology (IHuang et al.l. l2013l) and the ener ^ dissipation statistics to check the 
Lagrangian version rehned similarity hypothesis (IHuang and Schmittl . 120141 ). Here the 
inertial range 2 < r/y < 60 detected by MSA is due to the different scale dehnition. 


3.3. 2D turbulence velocity field 

Two-dimensional (2D) tu rbulence is an ideal model for the large scale movement of 


the ocean or atmosphere (iBoffetta and Eckd . 


2012 


; iKraichnan and Montgomery . 


1980 
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Figure 7. (Color online) Measured structure-function Sq{i) for the two-dimensional 
turbulent velocity. Due to the scale-mixing in SF analysis, there is no power-law 
behavior. 



Figure 8. (Color online) a) Measured qth-order 'Dq{t) from MSA. A dual-cascade 
power-law behavior is observed in the range 0.004 < £ < 0.04 for the forward enstrophy 
cascade, and 0.15 < £ < 1 for the inverse energy cascade, b) The curves compensated 
by the least square fitted scaling exponent ({q) to emphasize the power-law behavior. 
The vertical solid line indicates the forcing scale £f = 0.0614. 


Tabelind . 120021: iBouchet and Venaillel . 120121). We recall h ere briefly the main theoretical 


results of 2D turbulence advocated bv iKraichnanI fllOOTI ). 

The 2D Ekman-Navier-Stokes equation can be written as 

c^tu + u ■ Vu = —Vp + z/ u — au + fu, (12) 

in which u(x, t) = [m(x, t), u(x, t)] is the velocity vector, n is the fluid viscosity, a is the 
Ekman friction an d f„ is an exter nal source of energy inputting into the whole system 
at scale if = 1/kf (IBoffettal. 120071) . Parallelly, the vorticity a; = V x u equation is 

dtui + -|-u ■ Vo; = u ^ ~ + L- (13) 


To keep the whole system balance, two conservation laws emerge. The hrst one is the 
so-called energy conservation, inducing an inverse energy cascade from the forcing scale 
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Figure 9. (Color online) a) Fitted scaling exponents ^(g) for both the forward 
enstrophy (o) and inverse energy (□) cascades. The two solid lines represent 
respectively the Krainchan’s prediction ^(g) = g for the forward cascade and ^(g) = g/3 
for the inverse cascade, b) Measured singularity spectrum /(a). The inset shows 
the enlargement part in 0.2 < a < 0.5. A broader range of 1.17 < a < 1.45 and 
0.54 < f{a) < 1 for the forward enstrophy cascade indicates a stronger intermittency. 


if to large scales, i.e., i > if, which then leads to a five-third law in the Fourier space 
above the forcing scale, i.e., 

Eu{k) = , ka ^ k kf for the inverse energy cascade, (14) 

where Eu{k) is the Fourier power spectrum of u, is the energy dissipation by the 
Ekman friction, ka is the Ekman friction scale and kf = 1/if is the forcing scale. Below 
the forcing scale i < if, the enstrophy conservation law yields 

Eu{k) = C'rjl^^kr^, kf k kt, for the forward enstrophy cascade, (15) 


in which is the enstrophy dissipation by the viscosity z/ and ki, is the 
viscosity scale. This double-cascade 2D turbulence theory has been recognized as 
one of the most important res ults in turbulence since Kolmogorov’s 1941 work 
(jFalkovich and Sreenivas^ , 2006 ). 

For the last few decades, numerous experiments and numerical simulations 


have been devoted to veri 


(iBruneai 


Boffettal . 


and Kellavl. 


2005 


y the above mentioned fo r ward and inverse cascades 


Rutgersl. 


20071: iFalkovich and Lebedevl . 


1998 


2011 


Tan et ah 


Kellav et ah . 19981: [Bernard et al.l. 12006 


20141) with partially verification 


of the theory by Kraichnan. For example, iBoffetta and Mnsacchiol fl2010h performed 
a very high resolution numerical simulation, up to a grid number N = 32,768^. 
They stated that due to the scale separation problem, numerically the dual-cascade 
requires a very high resolution for verification. The inverse and forward cascades 
were o bserved for the third-order v elocity structure function as predicted by the 
theory (IBoffetta and Mnsacchiol . |2010L see Figure 3). For t he inverse cascade , it is found 
that the inverse energy cascade is almost nonintermittent flNa.m et ah 


2000 


Tan et al. 


2014j ). while for the forward enstrophy cascade, the intermittency effect still remains as a 


open question. This is because, according to the convergency condition, the structure- 
function requires a Fourier spectrum E{k) ~ k~^ with fi G (1,3), the /3-limitation 
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flHuang et al.l. 120lOh . Coincidentally, the Fourier scali ng exponent for the en s troph y 
cascade is fi > 3, see equa t ion ff^ and discussion in Ref. iBoffetta and Musacchiol fl2010l) . 
Theoretically, Nam et ah 1 2nnnh found that t he Ekman friction leads to an intermittent 
forward enstrophy cascade flBernarcy . 120001) . which has been verihed indire c tly b y 
studying the passive scalar held, instead of th e vorticity held (IBoffetta et ahl. 120021 ). 
More recently, this claim has been conhrmed bv lTan et al.l (120141) using Hilbert spectral 
analysis. A log-Poisson model without justice is proposed to h t the forward enstrophy 
cascade scaling exponent, see detail in Ref. iTan et al.l (l2014l ). However, to identify 
whether the forward enstrophy cascade is intermittent or not is still a challenge in the 
sense of data analysis. 

The D NS data for pr esent analysis is based on a fully resolved vorticity held 
simulation (IBohettal . 120071) . with an artihcially added friction coefficient. Numerical 
integral of equation (IT^ is performed by a pseudo-spectral, fully deali ased on a doubl y 
periodic square domain of size L = 2tt with N'^ = 8192^ grid points (IBohettal. 120071) . 
The key parameters adopted are z/ = 2 x 10“®, a = 0.025 and the energy input wave 
number fcf = 100 with very short correlation in time. The velocity held is be solved 
from the Poisson equation of the stream function -0, i.e., u = [dy'fi, —dx'fi]- Totally, hve 
snapshots with 8192^ x 5 = 3.36 x 1 0^ data points a re used for analysis. More details 
of this database can be found in Ref. IBohettal (120071) . 

The conventional SFs are shown in hgure[7] and no clear scaling range can be 
observed, neithe r any indic a tion o f the aforementioned two regimes. As discussed above 
and also in Ref. Tan et ah J2(114 ). such outcome can be ascribed to the scaling mixing 
in SF analysis (iHuang et al.l. 120101) . For comparison the results from MSA are shown 
in hgureEla), in which two regimes with diherent scaling relations appear, specihcally 
in the range 0.004 < i < 0.04 for the forward enstrophy cascade and 0.15 < i < 1 
for the inverse energy cascade. To emphasize the observed power-law behavior, the 
compensated curves are then displayed in hgureElb) by using the htted scaling exponents 
^(g). Two clear plateaus appear, conhrming the existence of the dual-cascade process 
in the 2D turbulence. The scaling exponents ^{q) are then estimated in these scaling 
ranges by a least-square-htting algorithm. Figure |9]a) shows the measured dual-cascade 
^(g). The theoretical predictions by equations (IT4l) (i.e. g/3) and (IT5l) (i.e. g) are 
indicated by solid lines. Note that the measured forward cascade curve is larger than the 
theoret ical one. A sim ilar observation for the Fourier power spectrum has been reported 
in Ref. iBoffettal (120071 ). which is considered as an influence of the fluid viscosity u. 

To detect the multifractality, the singularity spectrum f{a) is then estimated, which 
is shown in £gure|9]b). Based on the fBm case test, one can conclude that the inverse 
cascade is nonintermittent as expected; while the forward enstrophy cascade shows a 
clear sign of multifractality: a broad change of a and f{a). Note that a is the generalized 
Hurst number. The value range a G (1.17,1.45) implies a Hurst number H > 1, much 
larger than the one indicated by the equation (IT^. i.e., H = 1 . Suc h outcome may be an 
eff ect of the l o garith mic correction (IPasquero and Falkovichl . 120021) . while it is reported 
by iTan et al.l (120141) that the logarithmic correlation for the vorticity held is weak. It 
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also has to point out here that t he measured a and f (q) could be a f unctio n of z/ or the 
Ekman friction f Boffetta . 2007 : Boffetta and Ecke . 2012 : Tan et ah . 2014). Systematic 


analysis of the 2D velocity held with different parameter is necessary in the future for 
deeper insights. 


4. Discussion and Conclusions 


In summary, we propose in this paper a multiscale statistical method, namely multi-level 
segment analysis, without employing any decomposition. The statistical properties of 
extremal points at multiple s level inherit the intrinsic multi-scale physics of complex 
systems. For each s, the corresponding extremal point series dehnes the so-called 
segment, which can be characterized by the separation distance and the function 
difference at adjacent extrema. As an important extension of the existing work, the MSA 
method introduced in this paper proves meaningful in complex data analysis. This new 
approach has been verihed by revealing the monofractal property of the synthesized fBm 
processes, while the Hurst number H is slightly overestimated especially when H < 1/3. 
When applied to two numerical turbulent datasets, i.e. the high-resolution Lagrangian 
2D turbulence, MSA shows interesting outcome. The conventional methodologies, such 
as SF, fail to detect the inertial scaling of the velocity held, which is now recognized 
partially as the d-lim itation, and partially as contamination by the energetic structures 


flHuang et al.l . 120131 ). Very differently, MSA detects successfully the clear scaling ranges 


for both cases. More precisely for 2D turbulence, the retrieve d multifractal prop erty of 
the forward enstrophy conhrms the theoretical prediction bv iNam et ahl (120001 ). More 
systematic study of such multifractality and parametric dependance (e.g. the huid 
viscosity) is important in future to provide a better understanding of 2D turbulence. 

Finally, we provide the following general remarks: 


1 In principle, MSA is generally applicable without special requirements on the data 
itself, such as periodicity, the Fourier spectrum slope fi ( steeper than 3), noise 
perturbation and unsmoothness structures. 

2 The length scale in the context of MSA is determined by the functional structure 
rather than being an independent input. At different s levels the segments are 
different and thus the scales as well, which conforms with the multi-scale physics. 

3 A serious dehciency of the conventional SF is the strong mixing of different 
correlation and scaling regimes due to sample averaging, i.e. the hltering (as 
infrared and ultraviolet) effect. To dehne the segment structure helps to annihilate 
such hltering and extract the possible scaling relations in the respective scale 
regimes, which has been proved from analyzing the Lagrangian and 2D turbulence 
data. 

4 Moreover irregular sampling (e.g. missed points) is a common problem for data 
collection as in the geophysical experiments, LDV (Laser Doppler Velocimetry) 
measurement, etc. However, typically methods as Fourier analysis and others 
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require uniform spacial data points. Such trouble is easy to overcome in MSA 
since only extremal instead of all the functional points are involved. 
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